Phonon-mediated tuning of instabilities in the Hubbard model at half-filling 
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We obtain the phase diagram of the half-filled two-dimensional Hubbard model on a square lattice 
in the presence of Einstein phonons. We find that the interplay between the instantaneous electron- 
electron repulsion and electron-phonon interaction leads to new phases. In particular, a d 2 .2_ !/ 2-wave 
superconducting phase emerges when both anisotropic phonons and repulsive Hubbard interaction 
are present. For large electron-phonon couplings, charge-density-wave and s-wave superconducting 
regions also appear in the phase diagram, and the widths of these regions are strongly dependent 
on the phonon frequency, indicating that retardation effects play an important role. Since at half- 
filling the Fermi surface is nested, spin-density-wave is recovered when the repulsive interaction 
dominates. We employ a functional multiscale renormalization-group method- that includes both 
electron-electron and electron-phonon interactions, and take retardation effects fully into account. 
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I. INTRODUCTION 

The renormalization-group (RG) approach to interact- 
ing fermions^ has recently been extended to include the 
study of interacting fermions coupled to bosonic modes, 
such as phonons 1 . Experimental evidence indicates that 
in many strongly correlated systems, such as organic con- 
ductors and superconductors 3 , cupratesi^, cobaltates&, 
and conducting polymers^, both electron-electron inter- 
actions and phonons may play an important role. Also, 
recent advances in the field of cold atoms have made 
possible the creation of fermion-boson mixtures on ar- 
tificial lattices. In these mixtures the fermions interact 
through instantaneous on-site repulsion, and when the 
bosonic atoms condense, there is an additional retarded 
attractive interaction mediated by the fluctuations of the 
bosonic condensate; a scenario similar to that of elec- 
tron and phonon interaction in solid state systems^. The 
physics of the interplay between the repulsive electron- 
electron (e-e) and attractive electron-phonon (e-ph) in- 
teraction is not well understood and fundamental ques- 
tions arise that include a full understanding of retar- 
dation effects and whether competition/cooperation be- 
tween these interactions can lead to new phases. 

One of the most extensively studied model is the two- 
dimensional repulsive Hubbard model*«2iiLia, which in 
the absence of phonons and at half-filling becomes an 
antifcrromagnet due to nesting in the Fermi surface (FS) 
which drives an s-wave spin density wave (sSDW) insta- 
bility. In the presence of isotropic phonons and in the 
strongly coupled regime, it has been shown that the s- 
wave charge density wave instability dominates over an- 
tiferromagnetism and coexists with s-wave superconduc- 
tivity (sSC)i£. When anisotropic phonons are present 
and full retardation is taken into account the phase di- 
agram in parameter space associated with this model 
is summarized in Fig. Q, which reproduces the usual 
sSDW and sCDW instabilities but additionally a A x 2_ y 2- 
wave superconducting phase (dSC). Therefore, in this 
work we show that this dSC phase which appears at half- 
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FIG. 1: (Color online) Phase diagram for Einstein phonons 
of frequency uje = 1.0. Four phases involving antiferro- 
magnetism (sSDW) (purple squares), charge density wave 
(sCDW) (green stars) and s-wave (sSC) (blue circles) and 
d-wave (dSC) (red rhombs) superconductivity compete in the 
vicinity where the average phononic strength A approaches 
the bare on-site repulsion wo = 0.5. The lines distinguishing 
the different domains are guides to the eye. 



filling, is a result of cooperation between the two interac- 
tions. We stress that the role of the phonons goes beyond 
simply creating an effective attractive force between the 
electrons. Retardation plays a crucial role in generating 
the new phases, and the size of each region depends on 
the phonon frequency. This model can also be used to 
study various other systems where the e-ph interaction is 
present such as the quasi-2D k-(BEDT-TTF)2-X materi- 
als, which are stoichiometric with a fixed density of one 
electron per BEDT-TTF dimer— , or as we mentioned in 
the beginning, for fermion-boson optical lattice mixtures 
(where there is experimental control on the number of 
fermions per lattice point) deep in the bosonic conden- 
sate phase. This study is by no means complete in its 
results since it is limited to the weak coupling regime 
but is intended to chart the vast phase space of these 
type of systems with some preliminary results that will 



help arrange the important underlying physics and con- 
tribute to a better understanding of the competition and 
cooperation between the interactions involved. 

This paper is organized as follows. In section II we 
introduce the theoretical model and the RG method of 
analysis we employ. In section III we provide our results 
and draw the phase diagrams associated with the differ- 
ent orders in the system. In section IV we discuss the 
summary and the basic physics our work has highlighted 
and finally in the appendix we provide more details for 
the RG flow of the couplings and susceptibilities for the 
reader that is interested on the theoretical details. 



II. THEORETICAL MODEL 

The approach we employ in this study is based on a 
general RG analysis of a system of electrons recently 
expanded to involve the coupling of the electrons with 
phonons as welli. We use a generic model of electrons 
on a Hubbard lattice at half-filling interacting through 
the repulsive Coulomb interaction and being isotropically 
and anisotropically coupled to dispersionless bosonic ex- 
citations (Einstein phonons). The Hamiltonian associ- 
ated with this type of system is 

k,CT ki,k2,k3 cr 

+ U E J2 ^ + Z! 3(<7)4+q, CT C k, CT (6 q + 6-q) (1) 

q q.k 

where c^ (ck,<r) is the creation (annihilation) opera- 
tor of an electron with momentum k and spin u, and 
6„ (6 q ) is the corresponding creation (annihilation) op- 
erator of a phonon with momentum q. Also, £k = 
— 2£(cos k x + cos k y ) is the non-interacting electron energy 
at half-filling, lue is the Einstein frequency, uq is the e-e 
on-site repulsion, while the e-ph coupling g(q) is taken to 
be momentum dependent. Momentum conservation (up 
to reciprocal lattice vectors) implies k4 = ki + k2 k3. 
Going to the path-integral formulation, the bosonic fields 
can be integrated out exactly^ to give an effective e-e in- 
teraction 

LU 2 

U(h,k 2 ,k 3 ) =u - A| kl _k 3 | 2 E rj, (2) 



where ki = (iwi,ki), and X q = 2g(q) 2 /u>e- The bare 
interaction therefore contains an instantaneous on-site 
Coulomb repulsion term and a retarded attractive part. 
This is the interaction used in the functional RG analysis 
in this work. 



A. RG flow of the couplings 

The general idea of RG theory is to integrate out self- 
consistently the degrees of freedom in the electronic sys- 



tem in successive steps (called RG time) , so that the elec- 
tronic interaction becomes renormalized with each RG 
time step until it diverges, which is a signature of an in- 
stability. The approximate character of this method lies 
in the self-consistent evaluation of the renormalization 
of the electronic interactions, where bare values are as- 
sumed to be smaller in strength than the band-width of 
the available electronic energies. The same approxima- 
tion applies for the electron-phonon coupling strength as 
well (weak coupling) which results for the parameters of 
Eq. (Q to obey u ,ujE 7 g(q) <A = it. We also remain 
within one-loop accuracy^ which is adequate to capture 
the essential physics, but we expand the analysis to in- 
clude dynamics (frequency dependence) in the electronic 
couplings. In previous wor k 16 i 17 i 18 the couplings involved 
in the RG study were frequency-independent which led 
to self-energy corrections at one-loop that were neglected 
in order to keep the total number of electrons fixed. In 
this work, the implicit frequency dependence generates 
an imaginary part in the self-energy S^(fc) which we cal- 
culate in order to use the full cut-off dependent electron 
propagator, given by 



Ci{p) 



I£p|-A £ ) 



-£p-E<p(p)' 



(3) 



where 9{x) is the Heaviside step function, Ai is the value 
of the RG cutoff at the RG time £ = ln(A /A^) and 
£ p = ln(Ao/|£ p |), where Ao is the initial cut-off cor- 
responding to half the bandwidth (4t) as we defined 
above. As we have already mentioned the RG flow equa- 
tions are evaluated at one-loop^ accuracy using the gen- 
eral Polchinski equation 15 applied for this specific lattice 
model 16 . The zero-temperature RG flow equations for 
the couplings and the self-energy can be written asii 6 ^ 

d£Ui(ki,k 2 ,k 3 )= / -^[Ci{p)Ci{q pp - p)] 



xU i Jk 1 ,k 2 ,p)U i Jk 3 ,k 4 ,p) + / —[Ct(p)C t (p + q' ph )] 



d£ 



x[/~, (p,fci,fc 4 )#|/ (p,k 3 ,k 2 )+ / —[c e (p)Ce(p + q ph )] 
p h p h J p at 

-2U~ t (ki,p,ka)Ui (k 4 ,p,k 2 ) + U~ t (p,h,k 3 ) 



LI? (k 4 ,p,k 2 

'-ph 



Up (ki,p,ka)Ui (p,ki, 



(4) 



diZi{ki)=J - [Clip)] [Ueik uP ,p) - 2Ueik 1)P , h)} , (5) 

where we have defined q pp = k 4 + k 2 , q p h = k\—k 3 , q' h = 
k\ — K4, i pp = min\\t p , Zq pp — P f, t p h = min{z pi £p+q ph \, 
and £' ph = min{£ p , £ p + q > }. We have also used the short- 
hand notation f = f °° 1^- f ,„ ?■> . As we see, the above 

Jp J— 00 2tt J (2tt) ' 

RG flow equations are in general non-local in RG param- 
eter time £. 




RG time 



FIG. 2: (Color online) RG flow of the static homogeneous 
susceptibilities for the different instabilities for uq = 0.5, 
Ao = 0.6, Xtt — 0.4 and loe = 1.0. Inset: Corresponding 
RG flow of the couplings with uj\ — u>2 = u>3 — 0. The color 
scheme (according to decreasing divergence strength of the 
susceptibilities) is thick red line for dSC, thick purple line for 
sSDW, thick dashed cyan line for dSDW, thick doted blue line 
for sSC, thin dashed black line for dCDW, and thin green line 
for sCDW. 
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FIG. 3: (Color online) Phase diagram (mo = 0.0 and uje = 
1.0). Without on-site repulsion, the dSC phase disappears 
completely even for large phonon anisotropy. Color scheme is 
identical to Fig. Q. 



of completeness we provide the exact RG flow equations 
for the couplings and the self-energy. 



For the square lattice at half-filling the FS contains 
van Hove points where the density of states is logarith- 
mically singular. It is a well-justified approximation^ to 
divide the FS into regions around the singular points 
Qi = (tt, 0) + k and Q 2 = (0, n) + k with |k| <C tt 
(called two-patch approximation), where the majority 
of electronic states are expected to reside. This re- 
sults in four type of e-e scattering processes defined 
as 5i = t/(Qi,Q 2 ,Q 2 ), gi = £/(Qi,Qi,Qi), 53 = 
£/(Qi,Qi,Q 2 ), 54 = Ef(Qi,Q 2 ,Qi), which in the fully 
retarded case become frequency-dependent gi(toi, w 2 , 0)3). 
Deformations of the FS have been shown to be stable 
with respect to corrections due to e-e interaction o 11 ^? 
and need not be of a concern in this approximation. Dur- 
ing the RG flow the electronic density is kept fixed which 
results in the chemical potential flow canceling the flow 
of the real part of the self-energy. Non-locality in the RG 
equations is lifted because momentum transfers can only 
be zero or (71", n) and since £ k = — £k+(7r,7r) au RG param- 
eters map back to £. Phononic anisotropy is introduced 
by distinguishing e-ph scattering processes that involve 
electrons from the same patch (<7 2 , 54) and those that 
scatter an electron from one patch to another (51, 53), 
and by assigning different coupling strengths (Ao, A*-) to 
them. The only place that phonons enter into the RG 
flow is through the following initial conditions 



g{ t3 °(uJi,U!2,u>3) = uq - K 3 1 7-^ — — -^a, (6) 



UJ 2 E + (Wl - L0 3 ) 2 



g 24 (wi,w 2 ,W3) = ""o - A - 



LO% + (Wl - L0 3 ) 2 



(7) 



The reader interested in the analytic details of the flow 
equations can consult the Appendix where for the sake 



B. RG flow of the susceptibilities 

As we mentioned in the introduction, the different in- 
stabilities associated with the Hubbard model at half- 
filling are superconductivity (sSC, dSC), antiferromag- 
netism (sSDW, dSDW) , and charge density wave (sCDW, 
dCDW) with the corresponding couplings^ (we suppress 
the implicit frequency dependence for clarity) 



u (sd)SC = 


= -2(g 2 ±g 3 ), 




(8) 


u (sd)SDW = 


-- 2(54 ±53), 




(9) 


M (sd)CDW = 


= -2(2 5l ± 53 - 


-54), 


(10) 



where the signs are chosen so that strong fluctuations in 
a channel produce a positive value for the correspond- 
ing coupling, irrespective of its attractive (SC,CDW) or 
repulsive (SDW) nature. Since the couplings are func- 
tions of frequency and the most divergent ones are not 
necessarily at zero frequency^, only a divergent suscep- 
tibility can determine which phase has dominant fluc- 
tuations. The corresponding order parameters for the 
different phases are 



A SC (£,< 
A SDW (£,< 
A CDW (£,( 



Z^ 



Z^ 



Z^ 



a* ff , k (r)*_ CT! _ k (r), (11) 

*o- J k(r)*_ . )k+ ( WiW )(r), (12) 
*a,k(r)* ff , k+(7r!7r )(r), (13) 



and involve the creation of particle-particle (p-p) and 
particle-hole (p-h) pairs at given angle 6 and energy £ 
set by k. The homogeneous frequency-dependent suscep- 
tibilities associated with these order parameters are cal- 
culated by extending the one-loop RG scheme of Zanchi 





VtW 




FIG. 4: (Color online) Phase diagram for the same on-site 
repulsion (uo — 0.5) as in Fig. Q, but for a smaller phonon 
frequency (uje — 0.1). Slower phonons suppress supercon- 
ductivity over the nesting-associated CDW and SDW phases. 
Color scheme is identical to Fig. Q. 



FIG. 5: (Color online) Evolution of 52(^1, —0^1,013) coupling 
for Ao = 0.6, A,r = 0.4, and uje — 1.0. The RG steps cho- 
sen are I = 0.4, 2.4, 3.4, 4.9 (top left, top right, bottom left, 
bottom right). The color scheme is from lower dark blue (at- 
tractive) to higher dark red (repulsive) values. 



and Schulz^. Their definition involves scattering pro- 
cesses between pairs at different angles 9 with energies 
corresponding to fast modes and are given by 

x'WiMn -7*|) = n?=i/«&e(l6l -WfoA) 

x(A s (^,e 1 ,T 1 )A s ^ 2 ,e 2 ,T 2 )}, (14) 

where J(£, 9) is the Jacobian for the k — > (£, 0) trans- 
formation, and 5 =(SC,SDW,CDW). In the two-patch 
approximation 9 can take only two values and each sus- 
ceptibility becomes a 2 x 2 matrix that we diagonalize 
to extract the symmetric (s-wave) and antisymmetric (d- 
wave) eigenvectors of each corresponding order. The in- 
terested reader is urged again to refer to the Appendix 
where we provide the full expressions for the susceptibil- 
ity flow equations of all the different orders calculated to 
one-loop accuracy. 



III. RESULTS 

In Fig. J5J we show the numerical solution of the RG 
flow equations for the couplings with u)\ — u) 2 = W3 = 
(inset) and the corresponding static (u> = 0) homoge- 
neous susceptibilities for the case of Ao = 0.6, A^ = 0.4, 
lue = 1.0, and uq = 0.5. All parameters are expressed 
in energy units of It (a quarter of the band width) and 
for the numerical implementation the frequencies are dis- 
cretized into a total of N = 41 divisions up to a maximum 
value of w max = 6.0. The particular choice of Ao > A^ 
enhances the attractive BCS type pairing processes asso- 
ciated with g 2 , while suppressing the dominant repulsive 
nesting channels of 53 (Eqs. (|6l7|l ). and tilts the balance 
between the usually dominant sSDW and subdominant 
dSC phase. When reversed (A w > A ), the attractive 
channels of g^ combined with g\ lead to a CDW instabil- 
ity. 



The general phase diagram for fixed values of uq = 
0.5 and oje = 1.0 is shown in Fig. JJ), where the e-ph 
coupling is parametrized by the mean value A = (Ao + 
X 1T )/2 and relative anisotropy r = Ao/(Arj + A,r). Previous 
work on the phononic effects in this type of system was 
limited along r = 0.5 (isotropic phonons) and A > no 
and found the nesting-related dominant CDW competing 
with sSC, while dSC was suppressed^. In our generic 
study, extending to all possible configurations of coupled 
e-ph systems, we find that close to X = uq and along 
r = 0.5 there are four competing phases including dSC. 
Deep in the repulsive region antiferromagnetism prevails 
as expected (the counterpart of sCDW), but for r > 0.5 
type of anisotropy (Ao > A^), there is a large region of 
dSC dominance (the case of Fig. J2J is a point in that 
region). Another candidate for this parameter space is 
dCDW (charge flux phase)2i but we find that while this 
channel does get renormalized significantly (Fig.[3J), the 
suppressed g\ and (73 couplings undermine its strength. 

It should be pointed out explicitly that the generic 
characteristics of the phase diagram in Fig. (JTJ are inde- 
pendent of the actual value of the repulsive uq. The same 
phase diagram is always expected near A = uq and de- 
tails such as adding a next-nearest neighbor hoping term 
in the Hamiltonian or doping away from half-filling will 
only enhance our findings associated with dSC since all 
nesting-related processes will then be additionally sup- 
pressed. In the absence of Hubbard on-site repulsion 
(uq = 0), not only sSDW does not occur, as expected, 
but the dSC phase also disappears completely as shown 
in Fig. ® for the uje = 1.0 case. Therefore, the dSC 
phase is not being driven solely by anisotropic phonons, 
but by the combined effect of anisotropic phonons, on- 
site repulsion and nesting of the FS. 

In order to demonstrate the importance of retardation, 
we show in Fig. J3J the corresponding phase diagram for 
slower phonons of uje = 0.1, keeping uq = 0.5 fixed. We 



see that as compared to Fig. (I), the SC regions are sup- 
pressed and pushed towards larger values of A. This is 
a clear indication that the internal dynamics of the cou- 
pling functions extend beyond renormalizing the effective 
e-e interaction to attractive values. There is a rich inter- 
nal frequency structure among the couplings as we show 
in Fig. (J5J), where we plot the evolution of the pairing 
channel of 32(^1, — ^1,^3) for the case referring to the 
dSC instability of Fig. J5J • At the beginning of the flow 
we see the attractive part of the coupling (limited around 
lu 3 = u>i), to be much weaker compared to the repulsive 
parts. As the flow progresses, and we go beyond the ini- 
tial conditions of Eqs. (|6I7|I . the repulsive parts of the 
coupling are suppressed in w-space and in actual values, 
while the attractive parts proliferate in both. Remark- 
ably enough, when the critical point is reached we see 
that the same type of scattering process can become at- 
tractive or repulsive, depending on the dynamics of the 
electrons being scattered, with a range of strength that 
extends as much as three times the bare uq. This is a 
strong indication that retardation effects play a crucial 
and dominant role in these types of systems. 



such a system. In other words, what the functional RG 
study can beautifully provide to us is an unbiased study 
of the instabilities associated with this type of system in 
a self-consistent manner to the given one-loop accuracy 
order, which allows us the flexibility of probing an im- 
mense parameter phase space in a rather inexpensive way 
(numerically and analytically) . The assumptions of this 
approach are the weak coupling regime in the electron- 
phonon interaction and that the largest energy scale of 
the system is the bandwidth. Within this approximate 
method we are able to draw the first charts of the rich 
phase space when the fermion-boson interaction is no 
longer neglected but self-consistently included as well. 
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IV. CONCLUSIONS 

In conclusion, we have presented a one-loop functional 
renormalization group analysis of the Hubbard model 
at half-filling, self-consistently including the electron- 
phonon interaction. We find that different values of the 
e-ph coupling energy and anisotropy can tune the sys- 
tem into different instabilities, including sSDW, sCDW, 
sSC, and dSC. When only the instantaneous on-site re- 
pulsion is present, sSDW is well known to dominate. In 
the absence of on-site repulsion, phonon mediated attrac- 
tion drives the system into sCDW and sSC phases. In 
the presence of both interactions, depending on the com- 
petition between them, these phases appear in different 
parts of the phase diagram. In addition, these interac- 
tions also cooperate to generate a dSC phase. This phase 
therefore only emerges because of the interplay between 
the physics of Coulomb interactions and phonons. Re- 
tardation effects play an important role in the onset of 
these phases and also determine the size of the different 
regions in the phase diagram. 

Our work is not to be considered as a complete study 
of the two-dimensional Hubbard model in the presence 
of the e-ph interaction but as a preliminary but self- 
consistent study that hints towards the right physics of 
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APPENDIX A: COUPLING AND SELF-ENERGY 
FLOW EQUATIONS 

The renormalization-group flow equations for the four 
couplings and the self-energy are directly derived from 
Eqs. illlofl when the two-patch approximation is used 
which maps all momenta transfers either to (0, 0) or 
{■n, 7r) converting the intrinsically non-local flow equa- 
tions to local. This procedure is very well presented in 
the work of Zanchi and Schulai^ for the reader interested 
in the full exposition of details. Here we only highlight 
the basic points along with presenting the final formulas 
for the sake of completeness. 

In the two-patch approximation the whole FS is di- 
vided into two patches (each patch has its redundant 
mirror image). The RG flow Eqs. (111511 involve a delta 
function constraining the electronic energies £k to be "on- 
shell" which can be above {+^1} or below (— A^) the FS. 
This constraint induces a one-to-one correspondence be- 
tween k and {6, |£k| = A^) which consequently is em- 
ployed to simplify the 2D-k integration into an azimuthal 
integration over 9, with the proper Jacobian J(0, £k) in- 
troduced. The integrand is ^-independent and one can 
define a general operator F v acting on any product of two 
coupling functions gi(u>i, u>,w^)gj{oJ2, uj,ui4,) according to 



F v P[3i9j](wi,w,w 3 )(w2,w,w 4 ) = — 



_Ai r'\n W-At) 



d0- 



£(«>-Ai 



(w - E?(w)) (n - EJ'(fi)) + i/AJ 



+ OC 



dugi (aii , u , w 3 )gj (w 2 , w , ^4 ) 



[(w - E?(w)) (fi - E?(fi)) + i/A|] 2 + Af [u - EJ'(w) - u(Q - EJ'(fi))] 2 ' 



(Al) 



where fl depends on the external frequencies and u> (which is integrated out), v = ± and effectively defines two types 
of F operators, while E" is the imaginary part of the self-energy The Jacobian at half-filling has the convenient 
property <f J(9, —Ag) — <f> J(8, Ag). The complete RG flow equations for the four couplings and the imaginary part of 
the self-energy can then be written as 



{ui,U2,W3) = F + [u!+u>i -uJsli - 2[.g 3 .g3 + 9i9i](u>i, w,w 3 )(w 4 , u,uj 2 ) + [3353 + 9igi](^i, w,w 3 )(w, w 4 , ui 2 ) 



dgi 
d£ 

+ [5353 + 545l] (w,Wl,W3)(W4)W)W2)J +-F+[^ + CJ 3 -Wl]( -2[fl^flf 3 +Srifl'l](W3,W,a;i)(W2,W ) W4) 
+ [5353 +5154] (W3,W,Wl)(w,W2,W4) + [5353 +545l]( W 7 W 3,Wi)(w 2 ,W,W 4 ) I, 



(A2) 



dg2 
d£ 



(Wl,W 2 ,W3) = F_[Wl +CJ 2 - W]( [ff353 +52fir2](Wl,W2)W)(w3,W4,w) + [#353 + 0202] (u> 2 , U>1, w)(a>4, W 3 , w) I (A3) 



%3 
0^ 



(lji,w 2 ,w 3 ) = i^ + [w + wi -w 3 ] - 2[.g 3 .gi + giga](wi, w,w 3 )(w4, w, w 2 ) + [5354 + 5i53](wi, w,u 3 )(lj,lJ4, w 2 ) 



+ [5351 + 5453] (W, Wl> W3) (W4, W, W2) ) + ^+[^ + ^3 -Wi]( -2[g 3 .gi +.g!g3](LJ3,W,Wi)(w 2 ,W,W4) 



[5354 +5153] (W3,W,Wl)(w,W2,W4) + [5351 +34.g3](w,W3,Wi)(w 2 ,W,W 4 ) 

F_[a>i +cj 2 - w]( [5 2 g 3 + g?,g 2 ]{ui, u) 2 ,u>)(u> 3 , u>4,w) + [5233 +539 , 2](u>2,u>i,a')(w4,W3,w) 

F + [w + a;2 - w 3 ] [5354 + 5453] (w,w 2 ,W3)(w,W4,k>i) + -F+[w + w 3 - W2][ff354 + 54ff3](w,W3,W2)(w,wi,W4), (A4) 



— (wi,W 2 ,W 3 ) = F + [w + W 2 -CJ 3 ][53g 3 +g454](w,Lj 2 ,W3)(w,W4,Wl) 
I- F + [u + u 3 - W 2 ][ff353 + 5454](^,^3,W 2 )(w,Wi,W 4 ), 



(A5) 



<9^ 



M 



A, 



"VW-^ 



\4i 

I ,lk 



(0,-Ai) 



' d _ {[gi+92\(ui,uj,uj) - 2[g 2 +g 4 ](uj 1 ,uj,uj 1 ))(u - E"H) 
(w-££(w)) 2 +A| 



The above equations are numerically solved for each RG 
step £ until any coupling for any frequency channel di- 
verges to values greater than 20t at which point we stop 
the algorithm and form all the couplings associated with 
the different instabilities given by Eqs H8I1U|) . For fre- 
quency independent interactions the above equations re- 
duce to the usual^ 



APPENDIX B: SUSCEPTIBILITY FLOW 
EQUATIONS 



dgi 
d£ 


= -25i(5i - 


dg 2 
d£ 


9 2 

= -52-53 


#53 

d£ 


= -253(51 - 


dgi 
d£ 


2 1 2 

= 5 3 +5 4 - 



(A7) 

(A8) 

(A9) 

(A10) 



Once the RG flow for the couplings is numerically 
solved and the divergence point towards strong coupling 
is reached we calculate the susceptibilities associated 
with the major orders in the system. As we mentioned 
in the text, a general susceptibility calculation involves 
Eq. 1|14|) for all different orders. By including (to one- 
loop) all RG vertex corrections^ we obtain the general 
homogeneous susceptibility flow equation which in the 
two-patch approximation reduces to a 2 x 2 matrix given 



by 



A /" 7r / 4 f+oo 

n JO J-oo 



Ai r*i\ 

Wo " 

ZllZn + Z12Z21 Z\\Z\2 + Z12Z22 

Z21Z11 + Z22Z21 Z22Z22 + Z21Z12 



(Bl) 



where S = (AF, SC, CDW) and we have used the short- 
hand notation z i:j z lm = z 5 (9 l ,9 j ;uj 1 ,uj)z 5 (9i,6 m ;uj,uj 1 ) 
and the definitions 



Wt F {ui,u) = Wf DW (u>i,u) = 

[(^-S"(u;)) 2 +A2][( W +a; 1 -S^( W +w 1 ))+A £ 2 ] 



W e bt -(uJi,uj) = 

Ag-(o;-Sy(a;))(a;i-a;-Sy(a;i-a;)) 

[(^-E"H) +A2] [(^-w-E^-w)) +A2] 



(B2) 



(B3) 



The RG flow equations for the vertex functions associated 
with the nested related phases are given by 

A /" 7r / 4 f+oo 

d e zf F (cj 1 ,oj2) = -^ ddJ(9,-At) du>Wf F (ui,u) 



(B4) 



1" 3 7o 



21133 + 21254 zng 4 + z 12 g3 
22133 + 22234 22134 + 22233/ ' 



and 



d t zf DW (cj 1 ,uj 2 )=^rdej(e, -a,)/ duwf DH twi, w) 



(B5) 



^JO 



2ll3a + 2123/3 2 ll3/3 + 2123a 
2213a + 2223/3 2213,3 + 2223a/ ' 



where z^gk 



= „&( 



f,bJi,uj)gk{u)2,ijJi + u>,u>) and we 



used the additional definitions 



g a {uJ2, u>i+cj, to) = 2g 3 (u> 1 +w, w 2 , u))-g 3 (o<J2, wi-k*;, w), (B6) 
3/3(^2, wi+w, w) = 23i(wi-K^, W2, w)-34(^2, wi-fw, w). (B7) 



For the SC-related vertex function we have 



d iZ f c (uji,0J2) = m dSJ{B,-ht)\ duWf c (w u w) 



(B8) 



21132 + 21233 21133 + Z1232 
22132 + 22233 22133 + 2:2232 



where Zijgk 



= ,SCi 



i,6jl u>i,u>)gk(u>, wi -w, lo 2 ). All sus- 



ceptibilities are zero at the initial RG step I = 0, and the 
vertex functions obey 



2fco(^i> 6*2; ^1,^2) = -S 8u 9 2 S(uji - w 2 ). (B9) 



Also, due to the fact that we are at half-filling the flow 
equations for the nested related phases are locate. 
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